Temperature dependent magnetic properties of FePt: effective spin Hamiltonian 
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A model of magnetic interactions in the ordered ferromagnetic FePt is proposed on the basis 
of first-principles calculations of non-collinear magnetic configurations and shown to be capable of 
explaining recent measurements of magnetic anisotropy energy (MAE). The site (Fe,Pt) resolved 
contributions to the MAE have been distinguished with small Fe easy-plane and large Pt easy- 
axis terms. This model has been tested against available experimental data on the temperature 
dependence of MAE showing scaling of uniaxial MAE (Ki(T)) with magnetization (M(T)) K\(T) ~ 
M(T)"' characterized by the unusual exponent of 7 = 2.1. It is shown that this unusual behavior 
of the FePt can be quantitatively explained within the proposed model and originates from an 
effective anisotropic exchange mediated by the induced Pt moment. The latter is expected to be 
a common feature of 3d-5d(4d) alloys having 5d/4d elements with large spin-orbit coupling and 
exchange enhanced Stoner susceptibility. 

PACS numbers: 75.30.Gw 75.50.Ss 71.15.Mb 71.15.Rf 



INTRODUCTION 

Since the phenomenon of super-paramagnetism has 
been identified as one of the major limits for the con- 
ventional magnetic recording significant research ef- 
fort has been invested in the development of materials 
with large magnetic anisotropy energy (MAE). Recent 
progress in the fabrication and characterization of gran- 
ular and nano-particulate FePt films Q puts even more 
emphasis on the understanding of the giant MAE of FePt 
and its temperature dependence. The latter property ap- 
pears to be of critical importance for the development of 
future high density magnetic recording systems in par- 
ticular for heat-assisted magnetic recording 0. 

A systematic understanding of the temperature de- 
pendence of the MAE in itinerant magnets remains a 
challenge and one of the long standing problems in the 
theory of magnetism. The proposed model deals with 
mixed localized and itinerant magnetic moments and 
thus bears general importance as large anisotropy is 
achieved by combining strongly magnetic elements with 
non-magnetic ones, where the latter have large spin orbit 
coupling. 

The chemically ordered LIq phase of FePt has large 
uniaxial MAE with the first order anisotropy constant 
K x w 10 8 erg/cc Q based on the simple angular varia- 
tion of MAE E ams ~ Ki sin 2 6. In the Ll phase the cu- 
bic symmetry is broken due to the stacking of alternate 
planes of the 3d element (Fe) and the 5d element (Pt) 
along the [001] direction. It is well established that in this 
naturally layered ferromagnet the large MAE is mainly 
due to the contribution from the 5d element having large 
spin-orbit (s-o) coupling while the 3d elementprovides 
the exchange splitting of the 5d sub-lattice 

Theoretical framework for the temperature depen- 
dence of the MAE was mainly developed in the 50- 
60's. These efforts led to the formulation of a general 



quantum statistical perturbation theory (QSPT) sum- 
marized by Callen and Callen This theory pro- 
vides a general approach for calculating the magnetic 
anisotropy free energy for the effective spin Hamiltonian 
H = H lso + H anls with a large isotropic H lso and a rela- 
tively small anisotropic part H anls . The anisotropic part 
is taken in the form H an ' ts = ^ fc2-£g(Sj), where Lg(Sj) 
is a normalized polynomial of 2 nd order in the case of 
uniaxial symmetry with Si denotinga unit vector (in the 
classical case) at the atomic site i p ■ This form implic- 
itly assumes well localized magnetic moments leading to 
a universal parametric relation between MAE (K\(T)) 
and the magnetization M(T) Q. The two-sublattice 
modification of the QSPT proposed by Kuz'min Q has 
been successfully applied to describe MAE of the local- 
ized 3d-4f magnets. Note, that all these theories pre- 
dict that in the low temperature region MAE scales as 
K\ (T) ~ M 3 (T) except the very recent work by Skomski 
et. al. where a mean-field two sub-lattice Hamiltonian 
for Llo CoPt led to a K X {T) ~ M 2 {T) dependence @. 

Recent experimental results demonstrate that the uni- 
axial MAE of epitaxial FePt films can be very accurately 
fitted to a Ki(T)~ M 21 (T) dependence in the low tem- 
perature range 0, Il0| . This observation clearly demon- 
strates that the contribution of the single-ion anisotropy 
(leading to the M 3 dependence) is practically missing. 
Thus, the Hamiltonian used in the QSPT theory does 
not necessarily reflect all the essential features of the 
magnetic interactions in Llo FePt and possibly also of 
the other 3d-4d/5d ordered alloys. 

In the following we present a model of magnetic inter- 
actions in FePt which is constructed and parameterized 
on the basis of first-principles calculations and is shown 
to be capable of explaining on the quantitative level re- 
cent measurements of the K\ ~ M 7 (T) dependence with 
non-integer exponent 7 = 2.1 0,^3- Thus we propose 
a microscopic explanation of this unusual behavior and 
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test our microscopic model of magnetic interactions. 

The leading contribution to the anisotropic part of the 
spin Hamiltonian is described as anisotropic exchange 
mediated by the induced Pt atomic spin moments. The 
thermodynamic behavior of this Hamiltonian is investi- 
gated within the mean-field approximation (MFA) and 
in the classical limit using both Langevin dynamics and 
Monte Carlo simulations. We find that proper treatment 
of the magnetic interactions mediated by the induced Pt 
moment yield Ki(M(T)) and M(T) dependences in a 
good quantitative agreement with experiment, including 
the value of T c . Thus the proposed atomic scale model 
describes correctly the most important static magnetic 
properties and thus opens the way for modeling even 
more complex dynamic switching properties [TH ] . 

EFFECTIVE SPIN HAMILTONIAN MODEL 
BASED ON FIRST-PRINCIPLES CLSDA 
CALCULATIONS 

Our analysis begins with an investigation of the 
isotropic part of the spin Hamiltonian. We start 
with the constrained local-spin-density-approximation 
(CLSDA) calculations 01 f° r a non-collinear arrange- 
ment of Fe and Pt atomic spin moments as summa- 
rized in Fig. [U The CLSDA method [3 allows to re- 
duce the many electron problem to a minimization of 
the Hohenberg-Kohn energy functional Ehk(p(^), 
of charge, p(r), and spin density, a(r), with an ad- 
ditional constraint term which in the case of a non- 
collinear magnetic configuration leads to a CLSDA func- 
tional EcLSDA{p{'i f )-,S{r),hf') with an additional La- 
grange multiplier hf- having the meaning of an internal 
magnetic field |l2j. This magnetic field is determined 
self-consistently according to the condition of the de- 
sired orientation of the atomic moment rhi at the site 
i. The effect of thermal fluctuations on the electronic 
and spin sub-systems can be conveniently modeled with 
the spin-spiral (SS) configurations representing various 
degrees of short range order which is found in 3d mag- 
nets for temperatures well above T c [l4|. The values of 
the Fe and Pt moments for these SS configurations are 
calculated in the local coordinate system, associated with 
the orientation of the local quantization axis at site i as 
Mi — f_ F [n™ p (e) — nf n (s)\ de, where the local density of 
electronic states is a diagonal matrix over the spin indices 
n^' rT [lil Hf|. To summarize the most important results 
obtained for various SS configurations we introduce the 
convenient variable h v = H v /H® where H v is the ex- 
change field at site v of the Pt sublattice normalized by 
its value in the FM state H®. 

In Fig^we present the spin moments and total ener- 
gies of the non-collinear magnetic configurations calcu- 
lated self-consistently within the CLSDA, using a gen- 
eralization of the electronic structure method to treat 



non-collinear magnetic order In FigCJi we present 

normalized values of the Fe and Pt spin moments cal- 
culated as a function of h. The dependence shown in 
Fig. clearly demonstrates a dramatic difference in the 
degree of localization for Fe and Pt spin moments, respec- 
tively. The Fe spin moment remains almost constant as 
a function of h (or angle ff) indicating its relatively lo- 
calized nature in terms of the response to the thermal 
fluctuations. On the other hand, the Pt moment varies 
linearly with h. This result raises the important question 
how to develop an adequate model to describe a system 
with mixed localized and de-localized magnetic degrees 
of freedom. 

The theory of magnetic interactions due to localized 
magnetic moments is well established. Hence an effec- 
tive spin Hamiltonian associated with the localized (Fe 
moments) degrees of freedom can be constructed in the 
form 

ff toc = -XJ^S i .S j -^4^ ) [Sf] 2 , (I) 

which relies on configuration independent effective ex- 
change interaction parameters Jij and an effective single- 
ion anisotropy kp*J . The Si are used to denote Fe sublat- 
tice spin moments which can be treated as unit vector in 
the classical limit. The form of the spin Hamiltonian Eq. 
^ is well justified by our CLSDA results which clearly 
indicate that the Stoner excitations associated with Fe 
moments have much higher energy than those of the Pt. 

On the other hand, the value of the induced Pt mag- 
netic moment varies between and a maximum value 
for the FM state. The CLSDA total energy calculations 
without s-o coupling (isotropic energy) presented in Fig. 
^) allow to clarify how to deal with Pt magnetic degrees 
of freedom. Indeed, as the scalar-relativistic calculations 
show, the total energy associated with these delocalized 
degrees of freedom E l £° loc follows very closely the relation 
which can be derived from the Stoner-model expression 
for the total energy [J^ 

EZoc = de[n^{e) + n d "(s)} e - ±JJW* « 7[m,] 2 (2) 

where I is the intra-atomic exchange interaction param- 
eter, M® is the Pt magnetic moment in the FM state and 
n up , n dn are spin resolved densities of states forming delo- 
calized Pt moments denoted as M v with I v = 1/2I V [M®] 2 
and m v = M v jM®. The dependence shown in Fig. QJi 
also demonstrates that both the value and orientation of 
the Pt moment are entirely due to the exchange field of 
the surrounding Fe moments, following very closely the 
relation = x^H^, where \ v is the local Pt suscepti- 
bility constant. As follows from the linear dependence in 
Fig. ^a,b) , m v can be very accurately described by the 
relation, 
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FIG. 1: Results of the constrained self-consistent LSDA 
calculations for ordered Llo FePt: a) Fe (circles) and Pt 
(squares) spin moments normalized by their values in the FM 
state as a function of normalized exchange field h; b) isotropic 
(squares) and anisotropic (circles) part of the total energy as 
a function [m„] . Dashed lines used for the linear fit. 



where the Pt sublattice Xu is found to have a weak mag- 
netic configuration dependence; the Ji„ are the effective 
exchange parameters defined as the CLSDA total energy 
variation SEclsda/ 5^>i8n\ u in the FM state, where Si is 
the Fe and the Pt sub-lattice moment. The tempera- 
ture dependence of x„ arising from the Fermi distribution 
smearing is weak and will be neglected in the following 
statistical simulations. 

The anisotropic part of the effective spin Hamiltonian 
is calculated within the CLSDA, with s-o interactions 
included self-consistently, and is presented in Fig. ^p. 
We find that it can be be very accurately approximated 
by a quadratic dependence on the m„ parameter. We 
should emphasize that owing to the delocalized nature of 
5d/4d elements, in the general case this dependence can- 
not be guessed prior to the rigorous calculations. How- 
ever, our result for FePt allows us to identify the form 
of the spin Hamiltonian associated with delocalized mag- 



netic degrees of freedom, Hd e 
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(4) 



The magnetic energy is partitioned into localized and de- 
localized contributions using the CLSDA approach al- 
lowing for a unified description of the electronic degrees 
of freedom within the one-electron approximation. In 



particular, the Fe k F °J and the Pt single-ion k P u t ' con 
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previous study |4( we find that the Fe contribution is 
negative while Pt gives rise to a large easy axis contri- 
bution. The Fe and Pt contributions to MAE have been 
calculated within the LSDA and then corrected accord- 
ing to the previous LSDA+U calculations 5]. We find 



k [ p} = 1.427 meV and kg* = -0.097 meV which corre- 
sponds to the macroscopic uniaxial anisotropy constant 
K X {T = 0) = 7.7 10 7 e rg/cc. 

Finally, with Eqs. (Ill.'il 1|) we can introduce an effec- 
tive spin Hamiltonian reflecting all the above features 
revealed by our first-principles calculations. It is con- 
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structed as H = Hi 
convenient form 
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and can be reduced to the 
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We note that the spin Hamiltonian is now expressed 
in terms of the Fe degrees of freedom, with effec- 



tive exchange interaction parameters Ji 



Ji 



Xu \2 



Ylu JivJjv and an effective single -ion anisotropy, 
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and a two-ion anisotropy contributions 
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jL(0)/XfM3\7 t T. 
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As can be seen from these expressions, the Pt in- 
duced spin moments result in additional isotropic and 
anisotropic contributions, both depending on the effec- 
tive exchange interaction parameters Ji V defined in the 
ferromagnetic ground state. We find that J iv are rel- 
atively strong and positive, resulting in isotropic and 
anisotropic exchange interactions both stabilizing the fer- 
romagnetic order in the [001] direction (see [IS])- Un- 
like the Ji U: the effective exchange interaction parame- 
ters between Fe moments appear to be sensitive to 
the lattice spacing and may change sign from positive to 
negative as a function of the chemical ordering |l9l| . Im- 
portantly, the additional anisotropic contribution (EqsQ 
which controls Ki(M(T)) does not depend on J^. 

In order to assess the relative magnitudes of the single- 
and two-ion terms (Eq. |SJ , consider for clarity the nearest 
neighbor (NN) interaction only with N being the num- 
ber of JVAT's. Then the magnetic anisotropy free energy 
(Fanis(T)) within the first-order thermodynamic theory 
[I] (justified by (H anis ) / (H iso ) w 0.1, see Fig. HJd ) have 
an effective single and two-ion contributions with latter 
involving a sum over nearest neighbors 

F anis (T) a (H ams ) T = d^f\T) + (N- 1)4 2) / 2 (T),(8) 

where f^T) = (S l S i } T and f 2 {T) = (S , i S , j > T are the 
- d< 0) ' <0) ' 



tributions can be distinguished. In agreement with a single-site and pair correlation funct 



ions. 



u(0) 

^Fe 
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k$/N and df ] w k$/N within the AW approxi- 
mation. Given the small magnitude of K§\, the ra- 
tio between the single and two-ion contributions [N — 
1)S- 2 ^ /S ^ « (N — 1). Then, given the dominance of the 
two-ion contribution, and since within the MFA type 
approximation f 2 (T) ~ M 2 (T), one can arrive at a qual- 
itative explanation of the observed MAE temperature de- 
pendence. Clearly, for more accurate evaluation of the ra- 
tio between two and single-ion contributions the distance 
dependence of the 3i V has to be taken into account. In 
the following we present calculations beyond MFA and 
NN approximations to provide a quantitative analysis of 
the proposed model in terms of its ability to explain the 
non-integer exponent of the Ki ~ M 2A dependence. 



THERMODYNAMIC PROPERTIES 

In the following, we use Langevin dynamics simula- 
tions described in 20] to investigate properties of the 
effective spin Hamiltonian Eq. [5] within the classical ap- 
proximation. The exchange interactions are long-ranged 
and are taken into account for a distance of up to 5 
atomic unit cells via fast-Fourier-transformation, within 
the same calculation as the dipolar interaction. We sim- 
ulate spherical nanoparticles with open boundary con- 
ditions and sizes up to 14464 moments, corresponding 
to diameters up to 9.2nm. The anisotropy constant de- 
fined as the free energy difference between magnetiza- 
tion oriented parallel or perpendicular to the easy axis 
which according to the first order perturbation theory 
is given by the internal anisotropic energy difference 
Kx(T) = E a (T, B = BS±) - E a (T, B = Be\\) for the 
external field B . Then, according to Eq. single-ion 
and two-ion contributions can be distinguished as sum- 
marized in Fig. One can see that the two-ion term 
is the dominant contribution which is nearly nine times 
larger. Furthermore the data indicate a Curie temper- 
ature close to the experimental value of 750K [jj flf1 |. 
In Fig. [2t> we present the calculated K(M) dependence 
along with its single and two-ion contributions. Within 
these calculations, which we stress go beyond MF classi- 
cal approximation, the two-ion term scales as M 2 (T) in 
a wide T range, while the single- ion term follows M 3 (T) 
scaling only at low temperatures. Since K\(T) has both 
contributions, the expression for the low T expansion 



K x {T)/Kx{Q) « a M 2 + (1 - a) M 3 



M 



3-a 



(9) 



contains the a coefficient originating from the normal- 
ized two-ion contribution and the second term originat- 
ing from the single-ion contribution. The parameters d[°^ 

(2) 

and gQ_- following from our first-principles calculations 
allow us to evaluate finally the exponent of 3 — a = 2.09. 

Before proceeding to a comparison with experiments 
we examine the range of validity of our classical statis- 
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FIG. 2: (a) Ki(T) dependence using LD simulations with the 
effective spin Hamiltonian and its single and two-ion contribu- 
tions; (b-d) log- log plots for K i (T)/ K i (0) vs. reduced mag- 
netization M(T): (b) using LD simulation within the classi- 
cal approximation; c) calculated within the MFA QSPT for 
quantum j = 3/2 and classical j = oo, M 2 (T) and M 3 (T) 
dependence are presented for comparison; (d) comparision of 
K\(T) calculated with LD and experimental data of Thiele et 
al. and Okamoto et al. [Io||. The solid lines gives a power 
law with exponent 2.1. 



tical approximation. In Fig. we present results of the 
MFA QSPT calculations of the K\ (M) dependence and 
corresponding two and single-ion contributions. Both, 
the classical and the quantum K\{M) dependence are 
identical for the two-ion term in the whole range of tem- 
peratures and for the single-ion term in the low temper- 
ature range. Considering that the single-ion contribu- 
tion is dominant, we can compare our Langevin dynam- 
ics calculations with available experiment as shown in 
Fig. [2Ji- As one can see, our spin Hamiltonian with ab- 
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initio parameterization agrees very well in a wide range 
of temperatures, especially given that the low tempera- 
ture measurements by Okamoto, et al. also yielded an 
exponent of 2.1 [Io| . 

SUMMARY AND CONCLUSIONS 

To summarize, we propose an atomic-scale model of 
magnetic interactions in ordered Llo FePt with an ef- 
fective spin Hamiltonian constructed and parameterized 
on the basis of first-principles calculations. The pro- 
posed model is investigated analytically and using sta- 
tistical simulations. We find that the model describes 
on the quantitative level the experimentally observed 
anomaly in the temperature dependence of the magnetic 
anisotropy energy. We demonstrate that this observed, 
anomalous temperature dependence (Ki ~ M 21 (T)) is 
due to the delocalized induced Pt moments, leading to an 
exchange mediated two-ion anisotropy which dominates 
the usually expected M 3 contribution of the single-ion 
anisotropy. We believe that this mechanism is common 
for various 3d-5d/4d ordered alloys having 5d/4d nomi- 
nally non-magnetic elements with large s-o coupling and 
Stoner enhanced susceptibility. 

We thank R. Skomski, A. Shick, M. van Schilfgaarde, 
T. Schulthess, R. Sabirianov, C. Piatt and D. Weller for 
useful and stimulating discussions. 
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